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Preface 


Mathematics  has  and  will  always  be  a  fascinating  subject  for  me. 

The  almost  non-existence  of  methods  to  efficiently  reduce  matrices, 
especially  large  matrices  in  several  parameters  astounded  me,  since  a 
70-dimensional,  multi-parameter  state  matrix  is  commonplace.  My  main 
motivation  is  to  provide  a  tool  with  which  these  matrices  can  be  reduced 
to  a  more  sparse,  manageable  form,  while  maintaining  system  characteris¬ 
tics,  reducing  CP  time  and  memory,  thereby  reducing  cost.  This  thesis 
primarily  extends  the  work  of  Jones,  Leuthauser,  Frost  and  Storey,  ar.d 
Dolan. 
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Dr.  Dan  Repperger  of  the  AMRL  for  his  unfailing  computer  support  and 
expertise,  and  Professor  Albert  Moore  and  Major  Kenneth  Castor  for  taking 
time  to  review  this  report.  Most  of  all  I  would  like  to  thank  God  for 
his  never-ending  mental,  spiritual,  and  emotional  support  during  my  stay 
at  Wright-Patterson  AFB. 

Dedication  of  this  work  goes  to  my  thesis  advisor.  Dr.  John 
Jones,  Jr.,  whose  wealth  of  knowledge,  enthusiasm  and  moral  support  was 
indispensable  to  successful  completion  of  this  project.  Finally,  I 
thank  my  typist,  Jane  Manemann,  for  her  gargantuan  effort  in  bringing 
this  work  to  fruition. 


Charles  A.  Lew 
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Abstract 


^  Two  methods  of  reducing  state  matrix  equations  containing  several 
parameters  are  presented  in  this  thesis.  Pursuant  to  the  first  method, 
the  solutions,  obtained  by  iterative  and  non-iterative  methods,  of  the 
algebraic  Riccati  matrix  equation  and  Lyapunov  matrix  equation  are  relied 
upon  heavily.  The  second  method  hinges  on  an  equivalence  transformation 
to  a  Smith  form,  coupled  with  the  notion  of  zeros  of  a  matrix.  Main 
results  focus  on  the  two  parameter  case  but  can  be  extended  to  any 
number  of  parameters. 


v 


REDUCTIBILITY  OF  MATRIX  EQUATIONS  CONTAINING 
SEVERAL  PARAMETERS 

I.  Introduction 

Reducibility  of  the  equation 

x  31  Ax  (1.0) 


where  A  is  a  constant  matrix  has  been  widely  studied.  Classical  methods 
for  solving  Eq  (1.0)  involved  the  use  of  similarity  transforms  using 
eigenvectors  and  eigenvalues.  One  wishes  to  reduce  A  to  a  sparser  form, 
maybe  block  diagonal  or  even  diagonal  form  to  ease  computation  time, 
memory  requirements, and  cost  of  solving  Eq  (1.0).  Much  study  has  been 
done  for  the  case  where  matrix  contains  polynomial  elements  of  a  single 
parameter,  see  Jones  (Ref  13)  or  Browne  (Ref  8:138-149)  for  example. 
Presently,  there  is  no  known  method  for  reducing  matrix  equations  that 


have  polynomial  elements  in  many  parameters,  say  z^,z^, . . . ,z^,  or  z_  for 
short.  It  is  the  purpose  of  this  thesis  to  provide  a  method  for  reducing 


matrix  equations  in  several  parameters.  The  main  emphasis  will  deal  with 
the  two-parameter  case,  and  the  extension  to  the  multi-parameter  case 


can  be  done  by  induction  using  the  basic  results  found  here. 


Throughout  this  thesis  it  will  be  assumed  that  matrices  can  be 
rectangular  and  belong  to  the  ring  R[zJ.  Also,  conditions  of  stability 
and  nonsingularity  in  coefficient  matrices  will  not  be  Imposed. 

Two  methods  of  reduction  will  be  presented:  (1)  the  Riccati- 
Lyapunov  Method  and  (2)  the  Smith  Form  Method.  Iterative  and  non¬ 
iterative  solutions  of  the  Riccati  and  Lyapunov  equations  pertaining 
to  the  first  method  will  be  dealt  with  in  Chapter  II.  Chapter  III  shows 
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hOT  to  reduce  a  syste.  using  nechod  <»  -  (2),  and  gives  an 

application  to  singular  spate...  finally.  Chapter  IV  su_arites  what 
has  been  done  and  gives  recommendations  {or  future  study. 
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II .  Solutions  of  Rlccati  and  Lyapunov  Equations 

The  purpose  of  this  chapter  is  to  present  new  and  existing  methods 
for  obtaining  solutions  to  the  algebraic  Riccati  matrix  equation  and  the 
Lyapunov  matrix  equation.  Section  1  will  deal  with  the  Riccati  equation, 
presenting  four  methods  of  solution.  The  first  method  will  use  general¬ 
ized  inverses  in  its  solution.  The  second  method  will  implement  a 
function  of  a  matrix  polynomial  to  obtain  its  solution.  The  third 
method  draws  upon  a  power  series  expansion  to  arrive  at  its  solution. 
Finally,  the  last  method  utilizes  an  iterative  technique  to  zero  in  on  / 

a  solution.  Each  method  will  be  followed  by  an  example.  Section  2  deals 
with  solutions  of  the  Lyapunov  matrix  equation.  Since  the  form  of  the 
Lyapunov  equation  is  closely  tied  to  the  Riccati  equation,  the  methods 
developed  in  Section  1  will  also  apply  for  solving  the  Lyapunov  equation. 

An  additional  method  for  solving  the  Lyapunov  equation  that  will  be 
presented  is  the  Idempotent-Nilpotent  method.  Examples  will  follow, 
including  a  case  involving  one  parameter. 

Section  1.  Solutions  of  Riccati  Matrix  Equations 

The  Riccati  matrix  equation  occurs  most  commonly  in  modern  control 
theory,  in  the  design  of  optimal  controllers  (Ref  24:18)  and  estimators 
(Ref  29)  for  linear  time-invariant  systems  with  quadratic  cost.  Other 
uses  of  the  Riccati  equation  arise  from  singular  perturbation  methods 
(Ref  25)  and  filtering  theory  (Ref  20) .  New  uses  of  the  Riccati 
equation  crop  up  in  particle  beam  transport  processes  (Ref  6:28).  The 
form  of  the  steady  state  of  algebraic  Riccati  matrix  equation  we  will 
use  is  basically  the  following: 
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AX  +  XB  +  C  +  XDX  =  0 


(2.0) 


Other  equations  similar  in  structure  to  Eq  (2.0)  will  follow.  The  form 
of  Eq  (2.0)  actually  comes  from 

ATK  +  KA  -  KBR_1BTK  +  Q  =  0  (2.1) 

which  is  a  result  of  a  minimization  of  a  quadratic  cost  function  that 
is  related  to  the  linear  constant  dynamical  system 

x(t)  »  Ax (t)  +  Bu (t)  ,  x(tQ)  =  xQ  (2.2) 

Generalized  Inverse  Method .  The  emergence  of  generalized  inverses 
came  about  to  fill  the  void  in  the  solution  methodology  of  systems  with 
singular  or  rectangular  matrix  coefficients.  Applications  of  generalized 
inverses  include  iterative  methods  for  solving  nonlinear  equations, 
interval  linear  programming,  integral  solution  of  linear  equations 
(Ref  7:27,90,93,96),  estimation  theory  (Ref  10:1)  and  systems  science 
(Ref  21). 

For  every  finite  matrix  A  (not  necessarily  square)  whose  elements 
contain  coefficients  belonging  to  the  set  of  complex  numbers,  there 
exists  a  unique  matrix  X  satisfying  the  four  Penrose  equations  (Ref  7:7): 


AXA  - 

A 

(2.3) 

XAX  - 

X 

(2.4) 

(AX)*  - 

AX 

(2.5) 

(XA)*  - 

XA 

(2.6) 

where  A*  denotes  the  conjugate  transpose  of  A.  The  matrix  X  ■  A 
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satisfying  Eq  (2.3)  is  called  a  generalized  inverse  or  a  1-inverse, 


not  necessarily  unique.  The  matrix  X  =  A^,  not  necessarily  unique, 
satisfying  Eqs  (2.3),  (2.4)  is  called  a  reflexive  generalized  inverse 
or  a  1,2-inverse.  The  matrix  X  =  not  necessarily  unique,  satis¬ 
fying  Eqs  (2.3),  (2.4),  (2.5)  is  called  a  left-weak  generalized  inverse 
of  a  1,2,3-inverse.  The  matrix  x  =  A  ,  not  necessarily  unique,  satis¬ 
fying  Eqs  (2.3),  (2.4),  (2.6)  is  called  a  right  weak  generalized  inverse, 
or  a  1,2,4-inverse.  Finally,  the  matrix  X  **  A+  satisfying  Eqs  (2.3) 
through  (2.6)  is  called  a  Moore-Penrose  inverse.  The  matrix  A+  is 
unique,  and  is  the  same  as  A  \  the  ordinary  inverse  (Ref  7:7). 

Theorem  2 . 1 

X  is  a  solution  of  the  equation 

AX  +  XB  +  C  -  XDX  =*  0  (2.7) 


iff 


(A-I)K(LDK)rL  +  K(LDK)rLB  +  C  »  0  (2.8) 


where 


X  =  K(LDK)rL  (2.9) 

and  K  and  L  are  arbitrary  matrices  of  appropriate  dimension.  Equation 
(2.8)  is  called  a  consistency  condition  for  Eq  (2.7). 

Proof  (Necessity) 

Substituting  Eq  (2.9)  into  Eq  (2.7)  we  have 

AK(LDK)rL  +  K(LDK)rLB  +  C  -  K (LDK) rLDK(LDK) -  0  (2.10) 
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Using  Eq  (2.5)  we  have 


AK(LDK)rL  +  K(LDK)rLB  +  C  -  K(LDK) *  0  (2.11) 

(A-I)K(LDK)rL  +  K(LDK)rLB  +  C  =  0  (2.12) 


Proof  (Sufficiency) 

Expanding  Eq  (2.8)  we  have 


AK(LDK)L  -  K(LDK)L  +  K(LDK)LB  +  C  =  0  (2.13) 

AK(LDK)rL  +  K(LDK)rLB  +  C  -  K(LDK)  ^LDK  (LDK)  ^_L  =  0  (2.14) 

Comparing  Eq  (2.14)  with  Eq  (2.7)  we  have  X  =  K(LDK)rL. 

Method  for  Obtaining  A^ .  One  of  the  best  existence  theorems  of  a 
generalised  inverse  is  given  by  Jones  (Ref  14) .  The  following  treatment 
draws  upon  the  work  done  by  Jones.  Getting  the  reflexive  inverse  is 
easy  once  you  have  generated  two  distinct  generalized  inverses  or 
1- inverses.  Let  A^  and  be  two  1-inverses  to  find.  We  can  find  two 
non-singular  matrices  P  and  Q  such  that 


PAQ  = 


I  0 

r 

0  0 


(2.15) 


where  r  is  the  rank  of  A,  then  A^  is  given  by 


V  W 

VI  1 


(2.16) 


for  arbitrary  Ul,  VI,  W1  of  appropriate  dimension.  Similarly 
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p 

2  V  w 

2  2. 


(2.17) 


for  arbitrary  U2,  V2,  W2.  Now  that  we  have  two  1-inverses  we  can 
generate  A^  by  the  following  equation: 


A  *  A,  AA. 
r  1  2 


(2.18) 


See  Sontag  (Ref  27)  for  obtaining  the  generalized  inverse  of  matrices 
containing  parameters. 


Example  2.1 


2  -i  |  3  i]  r-i  i  i" 

B  »  C  -  D  -  (2.19) 

,-i  2j  Li  iJ  L-i  ->sj  li  l. 


Then  a  solution  to  AX  +  XB  +  C  -  XDX  *  0  is  given  by  X  -  K(LDK)rL  iff 
Eq  (2.6)  holds.  Choosing 


K  -  I  L  *» 


[H  h' 

L  *5  h. 


(2.20) 


we  have 


(A-I)K(LDK)rL  +  K(LDK)rLB  +  C  -  0 


(2.21) 


Hence 


X  =*  K(LDK) 


L  * 

r  U  h. 


(2.22) 
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As  an  extension  to  Eq  (2.5)  consider  the  generalized  Riccati 
matrix  equation  of  order  three  in  the  following  theorem: 

Theorem  2.2 

X  is  a  solution  of  the  equation 

AX  +  XB  +  C  -  XDX  -  XEXFX  =  0  (2.23) 

iff 

(A-I-K(LEK)  LF)K(LEK)  L  +  K(LEK)  L(B+DK(LEK)  L)  +  C  =  0  (2.24) 

r  r  r  r 

where 

X  *  K(LEK)rL  (2.25) 

Proof  follows  by  similar  methods  in  the  proof  of  Theorem  2.1. 

This  method  of  solution  can  be  extended  to  fourth  and  higher  order 
equations.  As  the  order  of  the  generalized  Riccati  matrix  equations 
increases,  so  increases  the  complexity  of  the  consistency  condition, 
placing  tighter  constraints  on  X.  This  implies  it  will  be  less  easy  to 
find  matrices  K  and  L  to  fit  the  consistency  conditions  corresponding 
to  their  appropriate  equations.  As  the  notion  is  rather  intuitive,  no 
proof  will  be  given. 

Example  2.2 


Let 

'  3/2 

1' 

'  -11/4 

-11/4' 

'1 

1' 

A  » 

.  1/2 

2. 

B  - 

.  “5/2 

-5/2  . 

C  - 

.  1 

1. 

'2  1' 

'  1  1* 

'  1 

l' 

D  - 

.  1  1. 

E  - 

.  1  1. 

F  - 

.  1 

3. 
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Then  a  solution  to  AX  +  XB  +  C  -  XDX  -  XEXFX  =  0  is  given  by 
X  =  K(LEK)rL  iff  Eq  (2.20)  holds.  Choosing 


K  =  I 


%  h 

L  hj 


(2.27) 


we  have 


(A-I-K(LEK)rLF)K(LEK)rL  +  K(LEK)^L (B+DK(LEK) ^L)  +  C  *  0  (2.28) 


Hence 


X  =  K(LEK)rL 


'  3s  ‘ 


(2.29) 


f (R)  Method .  This  method  relies  on  the  formation  of  a  state  matrix 
similar  in  form  to  the  one  in  the  Hamiltonian  system  of  equations. 
Referring  back  to  the  linear  system  in  Eq  (2.2),  one  can  form  the  control 
Hamiltonian  (Ref  2:238)  and  obtain  the  Hamiltonian  system  of  equations: 


X 

m 

r  -1  ti 

A  -BR  B 

T 

X 

-  H 

X 

y 

L -Q  -A  J 

Y 

Y 

L  J 

The  similar  system  we  will  use  is  the  following: 


$ 


D 

A  . 


$ 


(2.30) 


(2.31) 


Jones  (Ref  13)  has  an  excellent  paper  on  the  necessary  and  suffi¬ 
cient  conditions  concerning  solutions  to  Eq  (2.0).  Dolan  (Ref  11:65) 
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has  extended  Jones*  treatment  of  the  cases  involving  square  matrices  to 
the  rectangular  case. 


Theorem  2 . 3  (Jones  (Ref  13) 

Let  f  (X)  be  any  polynomial  of  degree  n  >  1  in  with  coefficients 
belonging  to  the  field  of  complex  numbers  such  that  for  square  matrices 
A,  B,  C ,  D , 


R  = 

'  -B 

D  ' 

II 

u/ 

■  u 

V  ' 

.  -c 

A. 

1  y  1 

.  M 

N  _ 

Then  a  solution  of 


(X,I)fy(R) 


with  U 


or  M 


-1 


existing,  or  a  solution  of 


(2.32) 


(2.33) 


f  (R) 
Y 


(2.34) 


with  M  ^  or  N  ^  existing  is  also  a  solution  of  Eq  (2.0). 

Experience  in  using  this  method  has  shown  that  Jones'  theorem  will 
not  always  work  even  if  U  ^  or  M  ^  or  N  ^  exists.  An  example  of  the 
inadequacy  of  Jones'  theorem  will  follow  after  slight  modification. 


Theorem  2.4 


Let  f^.(X)  be  any  polynomial  of  degree  n  >  1  in  \  with  coefficients 
belonging  to  the  field  of  complex  numbers  such  that  R  and  f(R)  are  as 
given  in  Eq  (2.32).  Then  a  common  solution  to  Eq  (2.33)  with  U  ^  or  M 
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existing  and/or  to  Eq  (2.34)  with  M  ^  or  N  ^  existing  is  also  a  solution 
of  Eq  (2.0). 

This  method  finds  all  solutions  that  Martensson's  method  (Ref  24) 
can  find  plus  more.  Illustration  of  the  theorem  is  followed  by  an 
example  from  Martensson  rewritten  to  conform  with  Eq  (2.0). 


Example  2.3  (Ref  24:31-32) 
Let 


o 

C\1 

'  2  O' 

O 

O 

1 

H* 

1 

A  = 

B  - 

c  = 

D  * 

o 

1 

o 

1 

.  0  3. 

L-i  -iJ 

(2.35) 


Then 


-B  D 
-C  A 


-2  0  -1  -1 
l 

0  1-1-1 
0  0  2  0 
0-3  0-1 


(2.36) 


and 


det(R-Al)  -  (A2-4)2 


(2.37) 


The  polynomials  f^,(X)  are  found  by  constructing  a  multiplication 
table  consisting  of  the  distinct  factors  of  the  characteristic  poly¬ 
nomial  ,  i .  e . , 
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Table  1 


Multiplication  Factor  Table  of  Eigenvalues  of  R 


X 

X-2 

X+2 

X-2 

(X-2)2 

X2-4 

X+2 

X2-4 

(X+2)2 

Examining  the  upper  or  lower  triangular  portion  of  the  table  we  find 
that  there  are  three  distinct  matrix  polynomial  cases  to  consider: 


Case  1.  R2  -  41 
~4  3  0  3 

0  4-30 

0  0  4  0 

0  0  3  4 


Case  1 
Case  2 


Case  3 


f  (R)  »  R2  -  41 
f  <R)  *  R2  -  4R  +  I 
f  (R)  *  R2  +  4R  +  41 


10  0  0 


fo  3  0  3 


-  4 


0  1 
0  0 


0 

1 


0 

0 


0  0-30 

0  0  0  0 


0  0  0  1 


0  0  3  0 


(2.38) 


'  U  M' 
V  N . 


(2.39) 


'0  0* 

r° 

-1/31 

■  0 

0  ■ 

»  -NM  «  - 

* 

.3  0. 

Ll/3  0  J 

.  0 

1 . 

0 

1 

H* 

U> 

- » 

‘0  3 

'  0 

o' 

X  -  M  U  - 

L 1/3 

oj 

.0  0 

* 

.  0 

1. 

(2.40) 


(2.41) 


12 


Case  2.  R2  -  4R  +  41 


16  3  4  7 

0  4-30 
0  0  4  0 

,  0  0  3  4J 


U  M 
V  N 


(2.42) 


X  -  -vu 


’o  0 

4  -3 

0  0 

1/64 

SB 

0  12 

0  16 

J 

0  -3 

J 

-1 


X  =  -NM 


X  -  MU  =  1/9 


'  0  o" 

4  -7 

0  0 

1/9 

=  1/9 

0  12 

.-1  \ 

-24  79 

4  -7 

16  3 

64  -16 

II 

-1  4 

o 

_ 1 

-16  13 _ 

(2.43) 


(2.44) 


(2.45) 


Case  3.  R  +  4R  +  41 


0  3-4-1 

0  12  -7  -4 

0  0  16  0 
0-12  34 


U  M 
V  N 


X  -  -NM 


-1 


16  0* 

-4  1 

64  -16 

1/9 

-  1/9 

3  4 

7  -4 

-16  13 

•  m 

r— 1 

1 _ 

0  3 

o 

o 
1 _ 

-  1/9 

- 

7  -4 

0  12 

0  -3j 

(2.46) 


(2.47) 


(2.48) 


Examining  the  common  solutions  to  the  above  three  cases  we  find  three 
solutions  satisfy  Eq  (2.0): 
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i 


i 


(2.49) 


•  -1 

* 

•* 

0  0 

0  0 

64  -16 

j* 

tt 

X2  = 

X,  -  1/9 

0  1 

0  -3 

J 

-16  13 

L 

L  J 

„  - 

Note  that  because  of  the  noncoramonality  of  X  in  Eq  (2.44),  it  is  not  a 
solution  to  Eq  (2.0)  as  Jones  claims,  even  though  M  ^  exists.  Also  let 
it  be  known  that  Martensson's  method  only  gave  solution  X^. 

The  next  theorem  is  a  slight  modification  of  a  theorem  by  Dolan 
(Ref  11:65),  who  treats  the  rectangular  case.  An  example  follows. 


Theorem  2.5 

Let  f  (X)  be  any  polynomial  of  degree  n  k  1  in  X  such  that  R  and 

f  (R)  are  as  given  in  Eq  (2.32)  for  A  ,  B  ,  C  ,  D  .  Then  a 
Y  6  n  mxm  nxn  mxn  nxm 

common  solution  to  Eq  (2.33)  with  U  ^  existing  and/or  to  Eq  (2.34)  with 

N  1  existing  for  X  is  also  a  solution  to  Eq  (2.0). 
mxn 


Example  2.4 
Let 


3  1 

-2 

A  =  [-1]  B  ■ 

C  *  [0-1]  D  * 

0  2 

L  j 

0 

m  * 

(2.50) 


Then 


-3  -1  -2 

0-2  0 
0  1-1 


(2.51) 


and 


det(R-Xl)  -  (X+3)  (X+2)  (X+l) 


(2.52) 


Using  Theorem  2.5,  two  solutions  to  Eq  (2.0)  are 
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X1  -  [0  1],X2  -  [1  0] 


(2.53) 


Beavers  and  Denman  (Ref  4:340),  using  a  similar  partitioning  scheme 

to  our  matrix  R,  claim  that  if  there  are  2n  distinct  eigenvalues  of  R 

then  there  are  „  C  =  (2n)!/(n!n!)  possible  solutions  to  an  nxn  matrix 
2n  n 

Riccati  equation.  Then  they  show  an  example  (Ref  3:141-143)  for  n  =  2, 
constructing  six  solutions.  Let  it  be  known  that  only  two  of  the  six 
claimed  solutions  actually  satisfy  the  matrix  Riccati  equation,  and  that 
those  two  solutions  were  readily  found  by  the  f (R)  method.  A  better 
proposal  to  the  number  of  possible  solutions  is  the  following:  For  any 
matrix  R,  if  there  are  p  distinct  eigenvalues  of  R,  then  the  number  of 
possible  solutions  is  given  by  p(p+l)/2.  It  is  easy  to  derive  this  number 
by  examining  the  Multiplication  Factor  Table  of  the  eigenvalues  of  R. 

Power  Series  Method.  This  method  relies  on  a  power  series  solution 
to  solve  Eq  (2.54)  and  Eq  (2.55).  As  the  quadratic  has  been  studied  by 
Dolan  (Ref  11)  and  the  cubic  case  by  Leuthauser  (Ref  18) ,  the  extension 
to  the  fourth  order  case  will  be  presented.  Note  that  the  higher  order 
case  can  be  reduced  to  the  Riccati  or  Lyapunov  equation  by  setting  the 
coefficient  matrices  of  those  terms  to  zero.  Consequently,  this  method 
can  be  extended  to  cases  involving  fifth  order  terms  to  finally  nth  order 
terms.  The  Eqs  (2.54)  and  (2.55)  arise  from  the  transport  equation  of 
particles  involving  three  scattering  processes  (Ref  28:633).  What 
follows  is  an  extension  of  Leuthauser  (Ref  18:59-64). 

Theorem  2.6 

Consider  the  equations 

A  /\  A  A  AAA  A  A  A  A 

AX-XB-C  +  XFX  +  XDXEX  +  XLXGXHX  (2.54) 
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A.  A. 


/V  A  A  A  A  A 


AY  -  YB  »  C  +  YFY  +  YDYEY  +  YLYGYHY 


(2.55) 


Equation  (2.54)  has  a  solution  X,  and  Eq  (2.55)  has  a  solution  Y  if  the 
equation 


A  A  AAA  A  A  A  A 


(A-XI)X(X)  -  Y (X) (B-Xl)  =  C  +  XFX  +  XDXEX  +  XLXGXHX  (2.56) 

have  solutions  X(X)  and  Y(X)  where 


X(X) 

=  xo 

+  Xxx 

+  X2X2 

+  .  .  .  + 

xpx 

p 

(2. 

.57) 

Y  (X) 

-  Yo 

+  XY 

+  x2y2 

+  ...+ 

XPY 

P 

(2. 

.58) 

5iven 

by 

A 

X  - 

xo  + 

XlB  + 

x2b2  + 

. . .  +  X 

p 

BP 

(2. 

59) 

A 

Y  = 

Yo  + 

AYi  + 

a2y2  + 

...  +  AP 

Y 

P 

(2. 

.60) 

To  obtain  the  X^  and  Y^,  we  must  examine  Eq  (2.56).  Equating  coefficients 
of  like  powers  of  X  we  have  the  following  set  of  equations: 

A  A  AAA  A  A  A  A 

C  +  XFX  +  XDXEX  +  XLXGXHX 

0 


(2.61) 


^0 

-  V 

AX, 

-  x„ 

-  Y,B 

+ 

Y„ 

1 

0 

1 

0 

AX„ 

-  X, 

-  Y„B 

+ 

Y, 

2 

1 

2 

1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

AX  , 

-  X  - 

-  Y  ,B 

+ 

Y  „ 

P-1 

p-2 

P-1 

p-2 

AX 

-  x  , 

-  Y  B 

+ 

Y  , 

P 

P-1 

-  X 

p 

P 

+ 

P-1 

Y 

P 

-  0 
-  0 
-  0 

2 

Multiplying  the  equations  in  Eq  (2.61)  on  the  right  by  I,  B,  B  , 
gives 
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> 


A  A  A.  A  A 


^0 

-v 

=  c 

AXLB 

-v 

-V2 

+  V 

=  0 

ax2b2 

-  xib2 

-V3 

+  Y^2 

*  0 

C  +  XFX  +  XDXEX  +  XLXGXHX 


AX  BP_1 

-  x  -BP_1 

-  Y  . BP  + 

Y  -BP_1  =  0 

p-1 

p-2 

p-1 

p-2 

AX  BP 

-  X  ,BP 

-  Y  BP+1  + 

Y  -  BP  =  0 

P 

p-1 

-  X  BP+1 

P 

+ 

p-1 

Y  BP+1  =  0 

(2.62) 


Columnwise  addition  gives 

[AX.  +  AX,B  +  AX_B2  +  ...  +  AX  nBP_1  +  AX  Bp] 

0  1  2  p-i  p 

-  [X-B  +  X.B2  +  ...  +  X  -BP_1  +  X  .BP  +  X  BP+1] 

0  i  p-2  p-i  p 

A  A  A  /V  A  A  A  A 

=■  C  +  XFX  +  XDXEX  +  XLXGXHX 


(2.63) 


Factoring,  we  have 


A(X_  +  X.B  +  X-B2  +  ...  +  X  -BP  1  +  X  BP) 

0  1  2  p-1  p 

-  (X„  +  X,B  +  ...  +  X  -BP~2  +  X  , BP_1  +  X  BP)B 
0  1  p-2  p-1  p 


A  A  A  A  A  A  A  A  A 

C  +  XFX  +  XDXEX  +  XLXGXHX 


(2.64) 


Thus  we  see  that  X  is  a  solution  to  Eq  (2.54).  Similarly  we  can  derive 

as  a  solution  by  multiplying  the  equations  in  Eq  (2.61)  on  the  left  by 
2  D+l 

I,  A,  A  ,  ...»  A1^  and  adding  columnwise.  This  gives 


A[I0  +  AY,  +  a2Y2  +  ...  + 

-  [T0  *  AY,  +  A2Yj  *  ...  a  AP"1Y(i_1  ♦  a'y^IB 

A  AAA  A  A  A  A 

-  C  +  YFY  +  YDYEY  +  YLYGYHY 


(2.65) 
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insuring  that  Y  is  a  solution  to  Eq  (2.55). 


Example  2.5 
Let 


r 

1 

1  0 

-1  -1 

!  36  -66 

A  = 

B  * 

c  = 

o 

rH 

- i 

0  -2 

-22  - 

-3J 

1 - 

o 

H 

i 

1 

o 

H 

1 _ 

0  0 

D  = 

E  = 

F  - 

= 

l  l 

.  0 

0  0 

=  L  =  G  =  H 


in  the  equations 


A  A  AAA 


AX  -  XB  =  C  +  XFX  +  XDXEX  +  XLXGXHX 

/S  A  A  A  AAA  A  A  A  A 

AY  -  YB  =  C  +  YFY  +  YDYEY  +  YLYGYHY 


Then  for  p  *  3 


X(X) 


Y  (X) 


1  1 
0  1 
-5  -11 
-3  2 


+  X 


+  X 


2  0 
0  0 
4  4 

1  -1 


+  X" 


+  X‘ 


0  0 

0  1 

-2  -1 
-1  1 


+  X' 


(2.66) 


(2.67) 

(2.68) 


+  X 


1  0 
0  0 
1  0 
0  0 


(2.70) 


Hence 


X  -  XQ  +  XjB  +  X2B  +  X3B 


Y  -  Yq  +  AYl  +  A2Y2  +  A3Y3 


-2  -8 
0  5 

-2  -8 
0  5 


(2.71) 


(2.72) 


Note  that  we  could  choose  as  many  or  as  few  terms  we  wish  by  adjusting  p. 
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Iterative  Method.  This  method  for  solving  the  Riccati  was  origi¬ 
nally  investigated  by  Leuthauser  (Ref  18:56-58)  and  Dolan  (Ref  11:37-41+) 
using  ideas  from  Lancaster  (Ref  17).  A  program  was  written  in  Fortran 
using  double  precision  arithmetic  and  implemented  on  the  CDC  6600. 
Computer  Processor  (CP)  time  per  iteration  was  very  small,  averaging  a 
miniscule  two  thousandths  of  a  second. 

Theorem  2.7  (Leuthauser  (Ref  18:50-58)) 

If  f(z)  =  (z+a) (z-a)  ^  where  a  ^  0  and  a  is  real,  and  if 

f (A)  =  (aI-A)_1(aI+A)  *  U  (2.73) 

f (B)  =  (al+B) (al-B)-1  =  V  (2.74) 

where  a  is  chosen  such  that  (al-A)  1  and  (al-B)  1  exist,  then  a  solu¬ 
tion  of 

X  =  UXV  -  ~  (U-I) (C+XDX) (V-I)  (2.75) 

2a 

is  also  a  solution  of 

AX  +  XB  =  C  +  XDX  (2.76) 

We  can  adjust  parameter  "a"  to  speed  convergence.  See  Leuthauser  (Ref 
18:56-57  )  for  proof.  If  conditions  in  Theorem  2.7  hold,  then  we  can 
obtain  a  recursive  formula  for  X: 

X  .  -  UX  V  -  (U-IKC+X  DXKV-I)  (2.77) 

n+1  n  2a  n  n 

Reinvestigating  the  examples  from  Martensson  (Ref  24:31-32+),  we  find 
that  we  can  obtain  more  solutions  than  Martensson  could  find  and  obtain 
a  greater  accuracy  and  faster  convergence  for  those  examples  than  could 
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Dolan  (Ref  11:71).  Rewriting  Martensson's  example  (Ref  24:44)  to  fit 
Eq  (2.0)  we  have  the  following: 


Example  2.6 


0  ll 

f 

I — 1 

H 

1 

o 

o 

A  =  B  *  j 

i 

c  = 

D  = 

.  1  0  - 

l  -l . 

o 

1 _ 

(2.78) 


Martensson's  eigenvalue-eigenvector  method  shows  two  non-negative 
definite  solutions  X^  and  X^  where 


Xl  = 


*  3W2  l-b/2  * 

x„  = 

‘  i/2-l  -/2+1  ' 

l+i^2  1+/2  _ 

2 

_  -/2+1  /2-1  . 

(2.79a) 


Dolan,  using  Eq  (2.77)  and  single  precision  arithmetic,  arrived  at  two  . 
solutions  X2  and  X^,  where 


*  .414 

-.414* 

X,  = 

* -2.41  2.41* 

.  -.414 

.414  . 

3 

2.41  -2.41  . 

(2.79b) 


using  an  initial  matrix  of  zero  for  both  solutions.  Solution  X2  was 
arrived  at  in  4  iterations  using  a  ■  3,  while  X^  was  obtained  in  8 
iterations  with  a  =  -7.  Although  Dolan  found  solution  X^  that  Martensson 
could  not  find,  he  could  not  find  solution  X^.  The  author's  attempts  at 
this  example  were  very  successful.  Three  solutions  of  Eq  (2.76)  were 
found.  Double  precision  arithmetic  and  sixteen  place  accuracy  were 
used.  Results  follow: 


*4.414213562373095  2.414213562373095* 

.2.414213562373095  2.414213562373095  . 


(2.80a) 
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using 


using 


4  2 

L  2  2 


.133 


.414213562373095  -.414213562373095 

L -.414213562373095  .414213562373095 J 


-1  1 
1  -1J 


,  a  =■  .2977 


X3  = 


-2.414213562373095  2.414213562373095 

2.414213562373095  -2.414213562373095 

-1-/2  1+/2  ‘ 

1+/2  -1-/2  . 


(2.80b) 


(2.80c) 


using 


-1  1 

L  1  -1 


,  a  *  -.1 


Number  of  iterations  to  converge  were  six  for  X^,  seven  for  and  eight 
for  X^.  Note  that  this  attempt  produced  the  intersection  of  solutions 
from  Martensson  and  Dolan.  One  can  see  that  for  different  values  of  a 
and  identical  starting  points  we  can  arrive  at  different  solutions. 

Also,  the  same  solution  can  be  arrived  at  using  different  starting 
points  and  different  values  of  a. 
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Dolan  (Ref  11:40)  gives  some  rules  of  thumb  that  work  quite  well. 
They  are  as  follows: 

(1)  Choose  an  a  equal  to  the  minimum  of  the  elements  of  the 
matrices  A,  B,  and  C. 

(2)  Assure  that  U  and  V  exist. 

(3)  Determine  if  the  sequence  generated  by  Eq  (2.77)  converges. 

(4)  If  the  sequence  is  converging,  then  iterate  until  less 
than  a  specified  tolerance. 

(5)  If  the  sequence  diverges,  increment  a  by  a  specified 

amount . 

The  specified  amount  in  step  5  could  be  determined  by  the  elements 
of  the  matrices  A,  B,  and  C.  For  example,  if  the  range  of  the  elements 
was  from  .1  to  10  then  an  adequate  adjustment  would  be  +  1.  But  if  the 
range  of  the  values  of  the  elements  was  from  1  to  100  then  an  incremental 
value  of  i  10  was  adequate. 

The  thrust  after  obtaining  an  a  to  get  convergence,  was  to  find 

a  .  ,  i.e.,  the  minimum  value  of  a  to  get  the  smallest  number  of  itera- 
min 

tions  for  convergence  within  a  specified  tolerance.  Experience  shows 

that  if  one  plotted  the  number  of  iterations  of  convergence  versus  a, 

the  plot  would  roughly  look  like  a  plot  of  cosecant  x  versus  x.  Also 

that  a  .  lies  in  a  certain  interval  and  is  not  unique.  For  example, 
min 

a  .  lied  in  the  interval  [.1329,  .133251  for  obtaining  the  six-iteration 
min 

convergence  of  X^. 

Dolan  did  not  specify  his  specified  tolerance  or  convergence 
criteria.  The  only  criteria  used  here  was  to  take  the  difference 
between  successive  approximations  and  see  if  it  Is  less  than  a  certain 
epsilon,  which  was  IE- 16  for  our  example.  Note  that  it  is  possible  to 
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converge  to  a  false  solution  using  this  difference  criteria.  One  must 
always  check  to  see  if  the  result  satisfies  Eq  (2.76).  An  example  of 
this  follows  from  Martensson  (Ref  24:31-32). 


Example  2 . 7 


A  =  B 


2  0 


0  0 


D  * 


-1  -1 


— 1 

rH 

1 

o 

LO 

-3 

I 

L-i 

Applying  the  iterative  method,  we 

find 

tht 

ec  results: 

o 

o 

■ 

"0  0 1 

X-  - 

X, 

3 

rH 

O 

- 1 

rH 

2 

. 

1 

o 

1 

u> 

i _ 

(2.81) 


(2.82) 


2.777777777777777  -1. 111111111111111 ' 

_ -1.111111111111111  1.444444444444444  . 


(2.83) 


Checking  Eq  (2.76)  we  find  that  X^  is  not  a  solution.  On  the  other 
hand,  if  we  apply  the  f(R)  method,  we  do  find  that  an  X^  does  exist, 
namely 


X 


3 


1/9 


'  64  -16  ' 
.-16  13  . 


(2.84) 


So  this  procedure  is  not  without  flaw.  Incorrect  solutions  can  appear. 

To  explore  bounds,  on  the  solutions  of  the  Riccati  matrix  equation 
to  get  an  initial  starting  point.  Is  not  of  prime . interest  here.  A  good 
treatment  of  bounds  may  be  found  by  Kwon  and  Pearson  (Ref  16)  and 
Bellman  (Ref  5) . 
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Section  2 .  Solutions  of  Lyapunov  Matrix  Equations 

The  Lyapunov  matrix  equation  arises  most  commonly  in  the  stability 
analysis  of  linear  and  non-linear  systems.  Other  applications  include 
the  analysis  of  beam  gridworks  (Ref  23).  The  numerical  solution  of 
certain  boundary  value  problems  in  partial  differential  equations 
(Ref  9),  and  the  design  of  optimal  control  systems  with  fixed  control 
structure  (Ref  19).  The  form  of  the  Lyapunov  equation  we  will  use  is 
the  following: 


AX  -  XB  *  C  (2.85) 

We  could  easily  solve  Eq  (2.85)  by  the  techniques  in  the  preceding 
section  by  setting  the  matrices  associated  with  the  quadratic  or  higher 
order  terms  to  zero.  Consider  the  following  example  with  one  parameter 
z,  using  the  power  series  method. 

Example  2.8 
Let 

T  1  1+8  0 1 


'  1  l+z" 

■-1-Z  -1-z  1  “ 

A  * 

B  - 

0  1  -z 

C  - 

.0  1  - 

0  l+z  z . 

L  0  0  1  J 

in  Eq  (2.85).  Then  for  p  *  1 


r 1  -1 

n  ri  o  oi 

X(X,z)  - 

+  M  ! 

(2.87) 

L-i  o 

ij  Lo  i  oj 

r 1  -1 

1 1  r1  o  °i 

Y(X,z)  - 

+  M 

(2.88) 

L-i  o 

l-z  J  1  0  1  0  J 
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Hence 


X  -  XQ  +  XjB  - 


z  z  1 
-1  1  1-z. 


(2.89) 


Y 


Y.  +  AY. 


'  2  z  1 
_-l  1  1-z 


(2.90) 


The  next  method  gives  an  explicit  solution  of  Eq  (2.85)  by  using 
the  principal  idempotents  and  nilpotents  of  the  coefficient  matrices  A 
and  B. 


Idempotent-Nilpotent  Method.  This  method  was  conceived  using  ideas 
from  Browne  (Ref  8:173,184-187,250-251)  and  relies  upon  the  decomposi¬ 
tion  of  the  coefficient  matrices  into  their  principal  idempotent  and 
nilpotent  components.  Proof  of  the  method  can  be  found  in  Jones  and 
Lew  (Ref  15) . 

Suppose  A  is  mxm,  B  is  nxn,  C  is  mxn,  and  A,  B,  C  belong  to  the 
complex  vector  space  of  p  by  q  matrices.  Let  A,  B  have  the  following 
repr esenta  t ion : 

m^<m  n^<n 

A  -  l  a  E  +  N  ;  B  -  £  b  F  +  >L  (2.91) 

j-1  3  J  J  k-1 


where  the  a^  are  m^  non-zero  distinct  complex  eigenvalues  and  the  b^ 
are  n^  non-zero  distinct  complex  eigenvalues.  The  {E^}  and  (F^)  form 
complete  sets  of  orthogonal  principal  idempotents  and  the  (N^ }  and  (Mk} 
form  complete  sets  of  orthogonal  principal  nilpotents.  See  Browne  for 
method  of  obtaining  idempotents  and  nilpotents.  The  following  condi¬ 
tions  hold  for  A  and  B  respectively: 
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m^<m 

*»2-v  J*1vi,Ilik’#,J'k 

EjNj  -  N  -  N^;  E^  -  -  0,  j  *  k 

(2.92) 

n^n 

Fk2  ■  V  I,  Fk  ■ I:  Vt  ■  °- k  11  *• 

k=l 

Fk“k  -  «k  •  W  FA  ■  Vi.  ■  °-  k  ^ 1 

(2.93) 

s. 

Theorem  2.8 

Let  A,  B  have  the  representation  given  in  Eqs  (2.91),  (2.92) 

» 

(2.93).  A  sufficient  condition  for  Eq  (2.85)  to  have  a  solution, 

is 

that  whenever  we  have  a  pair  of  indices  j,k  such  that  a^  *  b^,  then 

YFk  -  “j^k 

(2.94) 

EjCFk  - 0 

(2.95) 

> 

m^m  n^n  N  CF  -  E.CM. 

1  l  i  r  -  bJ- -  ~  ° 

j-1  k-1  aj  k 

(2.96) 

The  solution  X  is  given  by 

V®  Vn  E  CF. 

x-  l  l 

j-1  k-1  aj  Dk 

(2.97) 

i 
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Example  2.9 


Let 


2-11 

3  3-2 

4  10 


r -2  -1 
0  -1 


c  = 


0  0 
1  1 
1  1 


(2.98) 


Then 


A  =*  1 


2/4  2/4  -2/4 

0  0  0 

2/4  -2/4  2/4 

-1/4  5/4  -1/4 

+ 

10/4  10/4  -10/4 

+  3 

1/4  -1/4  1/4 

1 - 

u> 

4> 

CO 

h* 

4> 

1 

_ 10/4  10/4  -10/4 

3/4  -3/4  3/4 

a^E^  +  +  a2E2»  *  1;  ■  3;  ■  0 


(2.99) 


B  =*  -2 


1  1 

Lo  o 


-l 


o  -l 

Lo  l 


b^  +  b2F2;  bL  *>  -2;  b£  *  -1;  Mx  -  -  0 


(2.100) 


Finally 


„  .  E1CF1  ,  E1CF2  f  E2CF1  f  E2CF2 
al"bl  Vb2  a2“bl  a2~b2 


0  0 
1/3  1/3 

1/3  1/3 


(2.101) 


is  a  solution  of  Eq  (2.85). 
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III.  Reducibllity  of  Matrix  Equations  Containing 
Several  Parameters  with  Applications 

The  purpose  of  this  chapter  is  to  implement  the  knowledge  gained 
in  the  preceding  chapter  to  effect  a  block  reduction  to  any  given  state 
matrix  containing  polynomial  elements  in  several  parameters.  The  two 
parameter  case  is  dealt  with  here,  and  the  extension  to  n  parameters  is 
straightforward.  Two  methods  will  be  presented.  First  the  Riccati- 
Lyapunov  method,  then  the  Smith  Form  method.  Finally  an  application  to 
system  science  will  be  examined. 


Section  1.  Riccati-Lyapunov  Method 

This  method  from  Dolan  (Ref  11)  relies  on  solutions  to  the  Riccati 
and  Lyapunov  equations.  Dolan's  method  is  extended  to  the  rectangular 
case  with  parameters  . »Zq. 

Theorem  3 . 1 

The  differential  system  x  *  Ax  in  the  form 


(3.0) 


where  A  ,  B  ,  C  ,  D  belong  to  the  ring  of  matrices 
hlxou  nxn  mxn  tuuu 

R[z]  ■  R[z, ,z. , . . . ,z  ]  is  kinematically  similar  to 
—  12  q 


y 


D 

XD+A 


y 


(3.1) 


by  the  transformation 
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-X 


:] 


(3.2) 


where  X  is  a  solution  to  the  Riccati  matrix  equation 


AX  +  XB  +  C  +  XDX  =  0 


(3.3) 


Theorem  3.2 


The  following  differential  system  similar  in  form  to  Eq  (3.1) 


A  C' 
0  B 


(3.4) 


where  A  ,  B  ,  C  belong  to  R[z],  is  kinematically  similar  to 
nxn  mxm  nxm  —  1 


’  A  01 
.  0  Bj 


(3.5) 


by  the  transformation 


P 


(3.6) 


where  X  is  a  solution  to  the  Lyapunov  matrix  equation 


AX  -  XB  -  C 


(3.7) 


As  one  can  see,  the  method  is  a  two-stage  block  reduction.  Repeated 
application  of  the  method  to  each  block  and  sub-block  will  finally  result 
in  a  pure  diagonal  form. 

The  following  is  an  example  of  a  complete  reduction  for  the  constant 

case. 
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Example  3 . 1 


Reduce 


-3  -1  -2  -B  D 

0-2  0  x(t)  = 


01-1  -C  A 


From  Example  2.4  a  solution  to  Eq  (2.0)  is  X  *  [0  1].  Using  Theorem 

3.1,  Eq  (3.8)  reduces  to 


-3  -1/3  |  -2 


y(t)  =0-2  0  y (t) 


0  0  -1 


A  i  C 


0  B 


(3.9) 


Now  applying  Theorem  3.2,  Eq  (3.9)  reduces  to 


* -3  -1/3  O' 

z(t)  =  0-20  z(t) 


0  0-1 


*  A  O' 
.0  B. 


(3.10) 


using  X 


as  a  solution  to  Eq  (3.7). 


To  finish  the  reduction  we  apply  Theorem  3.2  to 


zx(t) 

ZjCt) 


-3  -1/3  zx(t) 


0  -2  J  |z2(t) 


(3.11) 


Using  X  *  1/3  as  a  solution  to  Eq  (3.7)  we  have 


r:  :]•■ 


(3.12) 
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Hence  from  Eqs  (3.9)-(3.12) ,  Eq  (3.8)  is  reduced  to 


-3  0  0 

z(t)  *  0-20  z(t) 

0  0-1 


(3.13) 


which  can  be  easily  solved.  Next  follows  a  reduction  in  two  parameters 
s  and  z. 

Example  3.2 
Reduce 


s  0  1 

1  10 

0  sz+1  1 

1  1-1 

! 

0  0  z 

z-s  z-sz-1  -1 

s  0  0 

(  ) 

0  sz+1  1 

0  0  z 

(3.14) 


Applying  Theorem  3.2  with  a  solution  to  Eq  (3.7),  namely 

'  s  0  o" 

X  =  0  10 

110 

we  have 


(3.15) 
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For  now,  Eq  (3.16)  will  be  left  in  upper  triangular  form,  even  though 
the  system  could  easily  be  solved  as  is  using  back  substitution.  Next 
a  reduction  method  using  an  equivalence  transformation  to  a  Smith  form 
in  two  parameters  s  and  z  will  be  presented.  The  following  treatment, 
marking  a  great  milestone  to  reduction,  is  found  in  Frost  and  Storey 
(Ref  12)  . 


Section  2.  Smith  Form  Method 

The  results  here  involve  transformations  of  equivalence  between 
pxq  matrices  A(s,z)  and  B(s,z)  of  the  form 

B(s,z)  “  M(s,z)A(s,z)N(s,z)  (3.17) 

where  M(s,z)  and  N(s,z)  are  unimodular  matrices.  The  equivalence 
transformation  of  most  interest  is  where  B(s,z)  is  a  type  of  Smith  form 
of  A(s,z)  over  R[s,z].  The  Smith  Form  method  for  two  parameters  is 
much  like  the  method  for  the  one-parameter  case  but  requires  no  zeros 
of  a  matrix  over  R[s,z]  for  equivalence. 

Theorem  3.3 

A  pxq  matrix  A(s,z)  of  rank  t  over  R[s,z]  can  be  reduced  by 
elementary  transformations  to  the  Smith  form 
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where 


f  E(s,z)  O]  ,  p  <  q 


S(s , z)  =  <  r 


E(s,z) 


VL 


E(s,z) 
0 


,  p  =  q 


,  p  >  q 


(3.18) 


E(s,z)  *  diag[e^(s ,z) ] 


(3.19) 


The  diagonal  elements  of  E(s,z)  are  the  invariant  polynomials  over  R[s,z] 
of  A(s,z),  given  by 


d^(s,z) 

ei(s,2)  =  d^Cs.z)  ’  1  =  1,2 . c 


(3.20) 


where  =  1  and  the  determinantal  divisor  (i  *  l,2,...,t)  is  the 
greatest  common  denominator  of  all  the  ith-order  minors  of  A(s,z),  and 
where  e^(s,z)  is  a  divisor  of  ei+^(s,z)  for  all  i  <  t.  That  is,  there 
exists  untmodular  matrices  M(s,z)  and  N(s,z)  over  R[s,z]  such  that 

M(s,z)A(s ,z)N(s,z)  *  S(s,z)  (3.21) 


Theorem  3.4 

Two  matrices  A(s,z)  and  B(s,z)  are  equivalent  over  R[s,z]  if  there 
exists  unimodular  matrices  M(s,z)  and  N(s,z)  such  that  Eq  (3.17)  is  true. 

It  is  easy  to  show  that  matrices  are  equivalent  over  R[s]  if  and 
only  if  they  have  the  same  Smith  form,  but  it  does  not  mean  that  matrices 
over  R[s,z]  that  have  the  same  Smith  form  are  equivalent.  See  the 
following  example. 
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Example  3.3 


s 

0 

1 

s 

0 

0 

A(s,z)  = 

0 

sz+1 

1 

B(s,z) 

35 

0 

sz+1 

1 

0 

0 

z 

0 

0 

z 

form  for  A(s,z) 

and 

3(s 

z)  is  given  by 

1 

1 

0 

0 

S(s,z)  = 

Vs 

,z)  = 

SB 

(s,z)  = 

0 

1 

0 

0 

0 

sz(sz+l) 

(3.22) 


(3.23) 


B(s,z)  and  A(s,z)  have  the  same  Smith  form,  but  there  does  not 
exist  M(s,z)  and  N(s,z)  such  that  Eq  (3.17)  holds.  Referring  back  to 
Example  3.2,  to  complete  the  reduction  via  Smith  form,  replace  block 
forms  by  form  in  Eq  (3.23).  To  get  a  stronger  condition  on  equivalence 
via  the  Smith  form,  the  definition  of  zeros  of  a  matrix  are  needed. 

Zeros  of  a  Matrix.  On  removal  of  the  determinantal  divisor  d^s.z) 
from  all  the  ith-order  minors  of  a  matrix  A(s,z),  the  remaining  poly¬ 
nomials  may  all  be  simultaneously  zero  for  one  or  more  values  of  the 
pair  (s,z).  Such  a  value  of  (s,z)  will  be  denoted  an  ith-order  zero  of 
A(s,z) . 


Example  3.4 

A(s,z) 


B(s,z) 


s+z  0  z 
0  s+z  0 
0  0s 
s+z  0  1 

0  s+z  0 
0  0  s 


has  the  same  Smith  form  as  the  matrix 


,  but  A(s,z)  has  first  and  second-order  zero 
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(0,0),  whereas  B(s,z)  has  no  first-order  zeros  and  many  second-order 
zeros. 

It  is  important  to  see  that  matrices  over  R[s,z]  which  have  the 
same  Smith  form  over  R[s,z]  but  do  not  have  the  same  zeros  are  not 
equivalent  over  R[s,z].  The  equivalence  transformation  over  R[s,z] 
preserves  the  zeros  of  a  matrix  over  R[s,z].  This  brings  us  to  the 
following  theorem. 

Theorem  3.5 

Two  matrices  A(s,z)  and  B(s,z)  are  equivalent  over  R[s,z]  iff 
A(s,z)  and  B(s,z)  have  no  zeros  and  have  the  same  Smith  form. 

The  following  is  an  example  of  reduction  to  Smith  form  using 
Theorem  3.3. 

Example  3.5 
Reduce 


using 


x(t) 


rs+1  1+z (s+1)  0 

s  sz+1  -(s+l)(s+z) 
0  s(s+l)  s+1 


x(t) 


M(s,z) 


1 

-1 


Ls2(s+1)  -s(s+1)2  1 


N(s,z) 


--Z  -[1+z (s+1)]  — (s+1) (s+z) [ 1+z (s+1) ] 
1  s+1  (s+1) 2 (s+z) 

0  0  1 


(3.24) 


(3.25) 
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We  have,  using  Theorem  3.3, 


y(t) 


1  0 
0  1 
0  0 


0 

0 


y(t) 


s+l+s(s+l)^(s+z)  J 


(3.26) 


The  Smith  form  works  well  with  matrices  containing  parameters  but 
one  has  to  work  hard  to  get  the  M  and  N  matrices.  Using  Theorem  3.3  as 
a  basis,  one  can  easily  extend  the  Smith  form  to  the  n-parameter  case 
by  induction.  This  is  shown  by  the  following  theorems. 


Theorem  3.3A 

A  pxq  matrix  A(z^)  of  rank  t  over  R[z_]  can  be  reduced  by  elementary 
transformations  to  the  Smith  form 


where 


(3.18a) 


E(z)  -  diagfe^z)]  (3.19a) 

The  diagonal  elements  of  E(z)  are  the  invariant  polynomials  over 
P.[zJ  of  A(z) ,  given  by 


d  (z) 

e^(z_)  m  j  ■  ^z")  *  ^  *  1,2, ...  ,t  (3.20a) 


where  d^  =  1  and  the  determinantal  divisor  d^(i  "  l,2,...,t)  is  the 
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greatest  common  denominator  of  all  the  ith-order  minors  of  A(z) ,  and 
where  e^z)  divides  e1+1(z)  for  all  i  <  t.  That  is  there  exists  uni- 
modular  matrices  M(^)  and  N(z)  over  R[z_]  such  that 

M(z)A(z)N(z)  -  S(z)  (3.21a) 

Theorem  3 . 5A 

Two  matrices  A(z)  and  B(z)  are  equivalent  over  R[zJ  iff  A(z)  and 
B(z)  have  no  zeros  and  have  the  same  Smith  form. 

Section  3.  Application  of  Reducibility  to  Singular  Systems 

In  this  section  we  will  illustrate  the  applicability  of  previous 
theorems  by  use  of  an  example  in  systems  science.  The  following  treat¬ 
ment  extends  Lovass-Nagy  (Ref  21). 

Consider  the  linear  system 

D(z)x(t)  *  A(z)x(t)  +  B(z)u(t)  (3.27) 

where  x  ,,  A  ,  B  ,  D  and  u  .  belong  to  R[z.  ,z„ , . . .  ,z  ],  or  in 
nxl  mxn  mxp  mxn  pxl  12  q 

shorthand,  R[zJ .  Equation  (3.27)  could  be  considered  as  the  state 
equation  of  a  control  problem,  where  one  hopes  to  find  an  input  u(t) 
that  will  force  x(t)  to  some  prescribed  function  of  t.  The  following 
theorems  are  necessary  to  analyze  Eq  (3.27). 

Theorem  3 . 6 

The  equation 

AXF  -  B  (3.28) 


has  a  solution  if  and  only  if 
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AA  BF  F  ■  B 


(3.29) 


and  the  general  solution  is  given  by 

X  -  A~BF~  +  W  -  A~AWFF~  (3.30) 

where  W  is  an  arbitrary  nxl  matrix. 

Theorem  3.7 

The  matrix  equation 


Ax  *  b 


has  a  solution  if  and  only  if 

(I-AA~)b  -  0 


(3.31) 


(3.32) 


and  the  general  solution  is  given  by 

X  -  A_b  +  (I-A_A)h  (3.33) 

where  h  is  an  arbitrary  nxl  matrix. 

An  example  of  the  reduction  of  Eq  (3.27)  follows. 


Example  3.6 

Let  A(z) ,  B(£),  D(z)  over 
a. (z)  c (z) 

A  - 

0  a2  (z) 


R[z]  be  as  follows: 


(3.34) 


where  a^,  c  and  b  are  polynomial  functions  in  the  parameters 
z^,...,z^  with  coefficients  in  the  field  of  real  numbers.  The  system 
in  Eq  (3.27)  can  be  reduced  to  an  equivalent  system  by  the 
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transformation  x  *  Py,  where 


(3.35) 


and  X  is  a  solution  to  the  Lyapunov  equation 


a^z^X  -  Xa2(z)  =  c(z) 


(3.36) 


Using  the  transformation  we  have 


DPy  =  APy  +  Bu 


(3.37) 


Premultiplying  by  P  gives 


P_1DPy  -  P_1APy  +  P~1Bu 


(3.38) 


Expanding  we  have 


1+X  -X-X 


1  -X 


a^z)  0 
0  a2  (z) 


1+Xb(z) 


(3.39) 


Dy  ■  Ay  +  Bu 


(3.40) 


Using  this  decoupled  system  the  procedure  to  find  u(t)  is  much 
simplified.  For  simplification  let  X  s  0.  It  is  easy  to  see  that  a 
generalized  inverse  of  B  is 


B  -  [1  0] 


(3.41) 


and  that 


(I-BB  ) 


0  0 


-b(z)  1 


(3.42) 
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The  consistency  condition  of  Eq  (3.31)  implies  that  Eq  (3.39)  has  a 
solution  u  if  and  only  if 


(I-BB-) (Dy-Ay)  =  0 


or  expanding 


If  y^  and  y^  are  any  two  functions  satisfying  Eq  (3.43),  then 
lar  solution  of  u(t)  is  given  by 


u(t)  =  5" (Dy-Ay) 


or  expanding 


u  (t) 


[1 


0] 


yrai(-)yi 

Va2(^)y2 


yl-al(-)yl 


(3.43) 


0  (3.44) 

particu- 


(3.45) 


(3.46) 


(3.47) 
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IV.  Conclusions  and  Recommendations 


Reducibility  of  matrix  equations  containing  several  parameters 
have  been  studied.  Two  reduction  methods  used  were  the  Riccati-Lyapunov 
method  and  the  Smith  Form  method.  The  Riccati-Lyapunov  method  used 
iterative  and  non-iterative  schemes  to  find  solutions  to  the  Riccati 
and  Lyapunov  equations.  This  method  to  effect  a  reduction  works  quite 
well.  Follow-on  work  should  be  the  formulation  and  comparison  of  com¬ 
puter  algorithms  that  demonstrate  the  non-iterative  solution  schemes 
delineated  in  Chapter  II.  For  the  iterative  scheme,  the  regions  of 
convergence  of  parameter  a  should  be  examined  more  closely,  the  method 
should  be  extended  to  include  cubic  and  higher  order  terms,  and  ways  of 
obtaining  bounds  and  initial  starting  matrices  should  be  explored.  The 
Smith  Form  method  used  an  equivalence  transformation  and  is  an  alternate 
form  of  reduction  that  can  do  the  reduction  in  one  step  as  opposed  to 
two  in  the  previous  method.  One  cannot  say  that  the  Smith  Form  method 
is  easier  than  the  other  because  they  are  equally  difficult.  As  it  is 
a  lengthy  process  to  obtain  solutions  of  the  Riccati  and  Lyapunov 
equations,  it  is  also  laborious  to  find  matrices  M  and  N  to  effect  a 
reduction.  Furthermore  the  effort  is  multiplied  when  the  matrices 
contain  multiple  parameters.  Nevertheless  the  Smith  Form  method  works 
very  well.  Future  work  concerning  the  Smith  Form  method  should  include 
the  formulation  of  a  computer  algorithm  that  will  reduce  a  given  pxq 
matrix  containing  several  parameters  to  Smith  form. 
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